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In the present work we present a fractal model for the j3 relaxation in supercooled liquids. The 
macroscopic dynamics is obtained by superposition of relaxation of independent mesoscopic regions 
(cages). Scaling relations and exponents are assumed for the distribution of cage sizes and for the 
size dependent response of the independent cages. In this way we obtain some scaling relations 
for the average response, in the time and frequency domain. For a particular choice of the scaling 
exponents we obtain the scaling relations predicted by the MCT. The comparison with recent light 
scattering data reveals that the scaling exponents of the distribution of cage sizes are universal but 
the single cage relaxation depends on the detailed material structure. 
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I. INTRODUCTION 



The liquid-glass transition at the calorimetric temperature T g has been a subject of extensive experimental and 
theoretical study for many years. However, more recently the interest has been shifted to the range of supercooled 
liquid temperatures, which is extended from T g to the melting temperature .J 1 ™. Light-scattering experiments in this 
range of temperatures reveals two relaxation processes, a and (3 relaxationscTLD . 

The commonly accepted explanation for the /3 relaxation is that it is related to the molecular details of the glass 
former. However, the mode coupling theory (MCTp suggests that it may be a general dynamical feature of dense 
disordered systems rather than being specific to the glass under investigation. Moreover, the MCT theoretically 
demonstrates the existence of the two distinct relaxation processes. The slower process (a) is associated with the 
diffusive motion and the faster one {(3) is attributed to strongly anharmonic motions of a limited number of molecules 
trapped in cages formed by their neighbors. 

On the other hand, computer simulations suggest that the (3 relaxation represents cooperative localized motions of 
a small number of molecules within a cage formed by neighboring moleculesEj. On the time scale of the /3 relaxation 
, these cages appear frozen. For longer times collective rearrangements continue allowing the eventual relaxation of 
the cages, which represents the a relaxation. 

If one neglects the interaction between different cages then one may model the j3 relaxation as the superposition of a 
heterogeneous distributiop_p£independent relaxing regions, the cages. There are several models in the literature which 
deal with such a situationliirEJ. In all cases, it is assumed that each cluster relaxes exponentially with a size dependent 
relaxation time. The main simplification of these models is thus the assumption of an exponential relaxation of the 
independent regions. The slow relaxation dynamics appears when one superposes the contribution of each region over 
the distribution of cluster sizes. However, the cooperative nature of the dynamics inside each cluster leads, by itself, 
to a slow relaxation dynamics. In fact, there are several models which -tLo-jnot contain disorder and obtain the slow 
relaxation as a consequence of the cooperative nature of the dynamicsali3E!3. 

In the present work we present a fractal model to describe the main features of the (3 relaxation in supercooled 
liquids. We assume that the macroscopic dynamics is given by the superposition of the relaxation of independent 
mesoscopic regions (cages). Nevertheless, different from previous works, we do not assume that each mesoscopic 
region relaxes exponentially. In this way we obtain some scaling relations for the average response, in the time and 
frequency domain. For a particular choice of the scaling exponents we obtain the scaling relations predicted by the 
MCT. The comparison with recent light scattering data reveals that the scaling exponents of the distribution of cage 
sizes are universal but the single cage relaxation depends on the detailed material structure. 

The paper is organized as follows. In section [n] we introduce the fractal model for the f3 relaxation in supercooled 
liquids and derive some of its predictions. The scaling relations for the normalized relaxation function and the 
imaginary part of the susceptibility are obtained. Then, in section [II we discuss the main results of the model, and 
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compare them with other theoretical approaches and recent light-scattering data reported in the literature. Finally, 
the conclusions are given in section IV. 



II. FRACTAL MODEL 



Computer simulations suggest that the (3 relaxation represents cooperative localized motions of a small number 
of molecules within a cage formed by neighboring moleculesEj. On the time scale of the (3 relaxation, these cages 
appear frozen. For longer times collective rearrangements continue allowing the eventual relaxation of the cages, 
which represents the a relaxation. In the following we restrict our analysis to the mesoscopic time scale between the 
microscopic characteristic time To and the mean live time of the cages, i.e. we focus our attention in the (3 relaxation 
regime. 



A. Single cage relaxation 



Let s be the number of molecules in a cage and t s the characteristic time associated with the cooperative rear- 
rangements of the molecules in a cage of size s. The relaxation in large cages, involving many molecules, is relatively 
infrequent and is associated with a long time relaxation process. On the contrary, the relaxation in small cages, 
involving a few molecules, occurs more often and is associated with a short time relaxation process. We thus assume 
that the characteristic time scales with the cage size according to 

Ts = T S C , (1) 

where £ > is a scaling exponent. In the following we measure time in units of tq. 

t s is the only relevant time scale in the p relaxation regime. Hence, the normalized response function will satisfy 
the scaling law 

$(s,t)=f>f 1 (s-h). (2) 

where fi(x) is a scaling function. On the other hand, has to satisfy the normalization condition 

nOC 

/ dt$(s,t) = l. (3) 
Jo 



This constraint carry as a consequence that 



r)=-l I dx^-=l. (4) 



B. Average relaxation 

The cage size follows certain distribution P(s). Computer simulations suggest that the number of molecules involved 
in cooperative rearrangements is widely distributcdO'lIa. Let s c ~ e~ 1 / cr be a characteristic cage size below with the 
distribution is scale invariant. We can therefore assume that P (s) satisfies the scaling law 

P(s) = s-\f2{e l/ °s), (5) 

where 1 < t < 2 and a > are scaling exponents, and f%(x) is a scaling function with the asymptotic behaviors: 
f2(x C 1) ~ const, and f2(% » 1) « 1. e is a control parameter, which characterizes how close is the system from 
the scale invariant case: P(s) ~ s~ r . It can be related to the mean cluster size through the equation 

(s) ~ e -( 2 - T )/ CT . (6) 

Hence (s) — >■ oo is equivalent to e — ■+ 0. In supercooled liquids and amorphous solids there is no long range order and, 
therefore, we expect that the mean cluster size is always finite, i.e. e > 0. We are thus considering the relaxation of 
scale invariant mesoscopic clusters. 
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The average normalized response function is given by 

$(t) = / dsP{s)®{s,t). (7) 
Jo 

Using eqs. (||), ([!]) and (pi) one obtains 

$(t)=t- 1 - a / 3 (e 1/A *), (8) 

where 

o=(T-l)r 1 , A = a/C, (9) 

/>0O 

fa(x)=C 1 dyy a - l h{y)h[{x/y) l l<). (10) 



/3(x) is another scaling function, expressed as a function of the previously introduced scaling functions /i(x) and 
/2(x). If feix) is analytical in x — then, expanding this function around x = 0, from eq. (Eo) one obtains 



/sC^wA-BiVC, (11) 

where 

A = C 1 ] HO) r dxx a ^- (12) 
Jo % 

5 = ^1/2(0)1 /"dra-VC^M. (13) 



jo J - 

Here we assume that the function/ 1 (x) is such that the integrals in eqs. (|4|), ( p"2| ) and ( |l3| ) converge. Substituting 
this asymptotic expansion in eq. (||) it results that 

$(t) - e (l+«)/A[( £ l/A^-l-a + C ( £ l/A^-1+^ (W ) 

where 

6 = r 1 -a, (15) 
and C = B/A. Eq. jl4| ) gives the asymptotic behavior of the response function in the scale invariant region t <C e _1 ^ A . 

C. Susceptibility 

The imaginary part of the normalized susceptibility %// is related to the normalized response function through the 
equation 

X»=FT"[*(t)]H, (16) 

where FT" denotes the imaginary part of the Fourier transform. Substituting $(i) by its scaling law in eq. (||) one 
obtains 

X>)=u; /4(e- 1/A w), (17) 
where /4(x) is another scaling function, given by 

r 





/ 4 (x) = / dycos(y)y- a - 1 f 3 (y/x). (18) 
Jo 

Moreover, in the frequency range e 1 ^ «C to -C 1 from eqs. (|l4]) and ( |l6| ) it follows that 



where 



CJ \ a , / UJ ^ ~ b 

+ b 



(19) 



X £ in ~ e a/A ^ min ~ e 1/A (20) 
From eq. ( p^ ) one concludes that the susceptibility has a minimum Xmin a ^ the frequency w m j n . 
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III. DISCUSSION 



In the previous section we have obtained the scaling laws for the normalized relaxation function and the imaginary 
part of the susceptibility. Now, we proceed to discuss these results and to compare them with other theoretical 
approaches and experimental data. n 

The scaling laws in eqs. (0) and (|l7|) are very similar to those obtained by the MCTeI. In fact, if we take 

A = 2a, (21) 

then we obtain the scaling laws of the MCT. This equation constitutes a scaling relation for the scaling exponents a 
and A, which has to be imposed to obtain the results of the MCT. Our formulation thus result more general, while 
the MCT prediction is contained as a particular case. 

From eqs. (§ and @ one can obtain the exponents r and £, introduced in our fractal model, as a function of a 



and b, which are usually measured in light scattering experiments in the (3 regime. It results that 

C _ 1 , r -*±£. ,22, 
a + b a + b 

Moreover, if one assumes that (^2|) is valid then form eq. @ we obtain 

2 a 



a + b 



(23) 

We emphasize again that eq. (|22j), and therefore eq. (|23j), is valid if one assumes the MCT prediction. Then using 
some light scattering reports for different supercooled liquids we have calculated these scaling exponents. The results 
are shown in | and D. 

First, we notice that the exponent r is larger than the one obtained in standard percolation (SP) in three dimensions, 
Tgp = 1.15. One of the shortcomings of standard percolation theory is that it neglects frustration, which introduces 
cancellations and hence the possibility of weak correlation, even if the interactions are much smaller than the thermal 
energy. Thus, standard percolation models are too .simple to explain the behavior of glass forming liquids. This 
conclusion has been also pointed out for spin glasseso. The inclusion of frustrated effects in percolation models is 
relatively receniO. The frustrated percolation (FP) model suggests that in glass-forming liquids there is a percolation 
temperature Tpp above T g , where frustration effects starts to be manifested. This transition leads to highly correlated 
localized clusters. 

Thus, we can describe the cage formation in supercooled liquids using a frustrated percolation model. Unfortunately, 
simulations in three dimensions are not yet available. In any case, it is expected that the exponents r and a are larger 
in FP than in SP. Frustration effects inhibit the formation of long range structures and, therefore, the distribution of 
cluster sizes in FP should decay faster than in SP, i.e. rsp < rpp. Moreover, the temperature T c may be identified 
with Tpp and e = |T — Tpp\. 

The data displayed in tables | and [n] has been divided in different groups attending to the structure of the glass 
formers. In the first class we have the two simple molecular glass formers ortophenyl and salol. These materials have 
a very simple structure with symmetric ( Van der Waals like ) interactions among their molecules. They are fragile 
glass formers and the scaling exponents a and b in both materials are quite similar. Molecular glass formers thus 
show similar static and dynamic behavior. Frustration may arise when local arrangements of molecules kinetically 
prevent all the molecules from reaching ordered close-packed configurations. Frustration effects have in this case a 
kinetic origin. In the second group we find the glass formers glycerol and B2O3. Glycerol is a fragile glass former 
which has a hydrogen-bonded network structure. B2O3 is a strong glass former and has a covalent bonded network 
structure. Hence, these materials are different in relation with the temperature dependence of the static viscosity. 
However, the scaling exponents a and 6, which characterize the dynamical properties, are in both materials quite 
similar. The difference in the temperature dependence of the static viscosity can be associated to the strength of 
the bonds. The hydrogen bonds are weaker than the covalent ones and, therefore, the network structure in grycerol 
is more sensitive to temperature shift than in B2O3. The common factor in these glass formers is the existence of 
network structures. We thus expect that frustration is mainly determined by the competition of different interactions 
generated by the multiple random bond network configurations. The existence of constrained bond angles or lengths 
leads to some optimized network structure, where some bonds will not satisfy these constrains and become partially 
broken (frustrated). Frustration effects have, in this case, a structural origin. Finally, in the third and fourth group we 
found the ionic glass former CKN and polymers, respectively. The exponents a and b are in both groups very different 
and greater than those reported for the first two groups. The origin of the frustration effects in these materials is less 
clear since we do not understand their structure well enough. 
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However, in spite that a and b change from one group to another, the scaling exponents of the distribution of cages 
sizes r and a are very similar. We found r w 1.34 and a « 1.69 (an exception are the polymers PPG and PET which 
are discussed below). Thus we have found a universal behavior for the distribution of cages sizes. On the contrary, 
the scaling exponent £ results material dependent. This scaling exponent is related to the single cage relaxation and, 
therefore, one may conclude that the single cage relaxation depends on the detailed structure of the glass former. The 
universal behavior can only appear when the single cage relaxation is averaged over the distribution of cages sizes. 

In the case of the polymers PPG and PET the exponent a deviates from that observed in the others glass formers. 
However, the scaling exponent r is similar to the one observed in the other glass formers. This deviation may be 
associated with the fact that the relation between r and the light scattering exponents a and b were obtained using 
general considerations, while that for a were obtained assuming that the MCT is correct. Hence, one may conclude 
that the MCT cannot be applied to polymers. Off course, to validate this conclusion a more extensive experimental 
analysis is needed. 



IV. SUMMARY AND CONCLUSIONS 



We have presented a fractal approach to the (3 relaxation in supercooled liquids. The main ingredients of the model 
are the existence of scaling laws for the single cage normalized relaxation function and for the distribution of cage 
sizes. In this way we have obtained the scaling laws for the average response, which for a particular choice of scaling 
exponents lead to those obtained by the MCT. 

We have obtained some scaling relations between the scaling exponents introduced in the model and those obtained 
from the fit to light-scattering data. Using these scaling relations and the reports for some light-scattering experiments 
we have found the (3 relaxation is indeed a general dynamical feature of dense disordered systems, which does not 
depend on the specific glass former under investigation. 

The scaling exponents distribution of cage sizes appears to be universal. The material dependence of the exponents 
a and b obtained from light scattering data is a consequence of the non-universality of the single cage relaxation, which 
depends on the detailed structure of the glass former. Moreover, the scaling exponent r and a were found larger than 
the one obtained in standard percolation in three dimensions. This discrepancy was attributed to frustration effects, 
and the cage formation was associated to a frustrated percolation problem. 

Our approach reveals that the MCT scaling relations may not be valid for polymers, but a more extensive analysis 
of experimental data is needed. The present work may be used as a guide for further experimental investigations of 
the P relaxation in supercooled liquids and glasses, in order to check our results. 
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Glass former 


a 


b 


c 


T 


Etef. 


salol 


0.33 


0.64 


1.0 


1.3 




orthophenyl 


0.33 


0.65 


1.0 


1.3 




B 2 3 


0.31 


0.62 


1.1 


1.3 




glycerol 


0.32 


0.61 


1.1 


1.3 




Ca .4Ko.6(NO 3 )i.4 (CKN) 


0.27 


0.46 


1.4 


1.4 




poly(propylene glycol) (PPG) 


0.23 


0.36 


1.7 


1.4 




poly(ethylene terephthalate) (PET) 


0.27 


0.44 


1.4 


1.4 





TABLE I. The scaling exponents £ and r for different supercooled liquids, obtained using the values of a and b reported in 
the corresponding references. 



Glass former 


A 


a 


salol 


0.66 


0.68 


orthophenyl 


0.66 


0.67 


B 2 3 


0.62 


0.69 


glycerol 


0.64 


0.69 


Ca .4Ko.6(N0 3 )i.4 (CKN) 


0.54 


0.70 


poly(propylene glycol) (PPG) 


0.46 


0.77 


poly(ethylene terephthalate) (PET) 


0.54 


0.76 



TABLE II. The scaling exponents A and a for different supercooled liquids. These exponents were obtained assuming that 
the MCT predictions are valid. 
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